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Q ■ Abstract 
00 

Stimulated by a scholium in Newton's Principia we find some beautiful re- 

^ ' suits in classical mechanics which can be interpreted in terms of the orbits in 

0\ ' 

\ the field of a mass endowed with a gravomagnetic monopole. All the orbits 

(N 



lie on cones! When the cones are slit open and flattened the orbits are ex- 
actly the ellipses and hyperbolae that one would have obtained without the 
gravomagnetic monopole. 

The beauty and simplicity of these results has led us to explore the similar 
problems in Atomic Physics when the nuclei have an added Dirac magnetic 



X 

H 

I monopole. These problems have been explored by others and we sketch the 

derivations and give details of the predicted spectrum of monopolar hydrogen. 

Finally we return to gravomagnetic monopoles in general relativity. We 
explain why NUT space has a non-spherical metric although NUT space itself 
is the spherical space-time of a mass with a gravomagnetic monopole. We 
demonstrate that all geodesies in NUT space lie on cones and use this result 
to study the gravitational lensing by bodies with gravomagnetic monopoles. 

We remark that just as electromagnetism would have to be extended be- 
yond Maxwell's equations to allow for magnetic monopoles and their cur- 
rents so general relativity would have to be extended to allow torsion for 



general distributions of gravomagnetic monopoles and their currents. Of 
course if monopoles were never discovered then it would be a triumph for 
both Maxwellian Electromagnetism and General Relativity as they stand! 
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I. INTRODUCTION 



One of us was asked to review Chandrasekhar's (1995) book on Newton's Principia (1687) 
for Notes and Records of the Royal Society (Lynden-Bell 1996). This led to reading passages 
of Cajori's translation of Principia. In his first proposition Newton shows that motion under 
the influence of a central force will be in a plane and that equal areas will be swept by the 
radius vector in equal times. In his second proposition he shows that if a radius from a point 
S to a body sweeps out equal areas in equal times then the force is central. There follows 
this scholium: "A body may be urged by a centripetal force compounded of several forces; 
in which case the meaning of the proposition is that the force which results out of all tends 
to the point S. But if any force acts continually in the direction of lines perpendicular to 
the described surface, this force will make the body to deviate from the plane of its motion; 
but it will neither augment nor diminish the area of the described surface and is therefore 
to be neglected in the composition of forces." 

What does this mean? 

The words described surface have been translated from a Latin word that carries the 
extra connotation of a surface described by its edge. We shall take this to be the surface 
swept out by the radius vector to the body that is now describing the non-coplanar path. A 
force normal to this surface at the body must be perpendicular to r and v which are both 
within the surface, so Newton is considering extra forces of the form NniQr x v where N 
may depend on r, v, i etc. We write the equation of motion 



modh/dt^ = -V'{r)r + NL, 



(1.1) 



where V{r) is the potential for the central force, f is the unit radial vector, and 



L = moT X V . 



(1.2) 



Taking the cross product fx (1.1) we have 



dL/dt ^ NrxL 



(1.3) 
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from which it follows either geometrically a la Newton or by dotting with L that 

|L| = constant. (1.4) 

Now if ip is the angle measured within the described surface between a fixed half-line ending 
at S and the radius vector, 

irV = ||L|/mo (1.5) 

so equal areas are swept out in equal times just as Newton says. To see this angle more 
precisely it is perhaps worthwhile to work in axes which are continually tilting to keep up 
with the plane of the motion. In any axes rotating with angular velocity 0(t), the apparent 
acceleration r is related to the absolute acceleration d^r/dt^ by 

d'^r/dt^ = r + 2fixr + f2xr + rXx(12xr). (1.6) 

We shall apply this formula to axes which are always tilting about r such that in these axes 
the motion appears as planar. Thus putting 17 = Qf in Eq. (1.6) 

d^r/di^ = r + Qr^^L/mo . (1.7) 

Inserting this into Eq. (1.1) and choosing 

Vt = rN, (1.8) 

we recover in these axes the equation we would have had in inertial axes had Newton's extra 
force oc N been absent i.e., 

mor^-V'r. (1.9) 

Thus relative to these moving axes r x mor = L is constant not only in magnitude but also 
in direction and 

\r X r\ = r'^ip = L/mo (1-10) 

where (p is the angle at S between some line fixed in the moving axes and the current radial 
line (this is of course equal to the earlier angle since this moving plane is 'rolling' on the 
described surface about the common radius vector). 
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We now return to the inertial axes in which the direction of L varies in accord with 
(1.3). Dotting Eq. (1.1) with v = dr/dt the N term goes out so the energy equation is left 
unchanged and we have, remembering that is constant, 



2 2 



^2 _^ (—\ r"^ 
Vmo/ 



+ V = E. (1.111 



Here r is the same in fixed or rotating axes since this r is scalar. Equation (1.11) demon- 
strates that the radial motion r{t) is precisely that which would have occurred had N been 
zero. Furthermore (1.9) and (1.10) demonstrate that within the tilting axes, or [using (1.5)] 
within the described surface, the solution r{(p) is precisely the same function that we would 
have found for the truly planar motion that occurs with N absent. Although this extension 
of Newton's theorem is not in Principia it would surprise us if Newton had not seen and 
understood it. There is interesting historical research to be done here on Newton's surviv- 
ing manuscripts. We know from Whiteside that this scholium was not in the first draft of 
Newton's De Motu Corporum written in Autumn 1684 but appears in its revision which is 
probably dated to the Spring of 1685. 

Although (1.5) and (1.11) are sufficient for the solution of the motion within the described 
surface, we need to find that surface by solving (1.3) for a complete description of the motion. 
This is not particularly simple and to do it we need to prescribe how N depends on r, v, t 
etc. However dL/dt and dr/dt are always parallel since both are perpendicular to r and L. 
This led us to consider under what circumstances they might be proportional. In particular 

dv/dt = -f X (f X v)/r = -r X L/(r^mo) (1-12) 

so in full generality we have from (1.3) 

dL/dt = -{moNr^)dr/dt . (1.13) 

This demonstrates that when moNr^ — — constant we have a beautifully simple solution 
to (1.13) to wit. 

L (5*r = j = const (1-14) 
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here j is the vector constant of integration; notice that has the same dimensions as L. 
Since L and r are perpendicular we deduce, dotting with r 

j-f = g, (1.15) 

which shows that the angle between j and f is constant so f moves on a cone whose axis is 
along j. Similarly dotting (1.14) with L we find = j • L so likewise L moves on another 
cone with j as its axis. If this cone has semi-angle x then L/|j| = cosx, but f and L are 
orthogonal and by (1.14) they are coplanar with j so we may choose cf. (1.15) 

a/|j|=sinx (1.16) 

so the angle between f and j is 7r/2 — x as shown in Fig. 1. Notice from (1.16) that the 
angle of the cone is determined completely from |L| and the force constant Q*. Orbits 
with larger |j| have smaller x so the angular momentum then moves around a narrow cone 
and r then moves around a very open one. For |j| ^ it is almost planar. Fig. 1 
illustrates two circular orbits moving in opposite senses about the same axis. Notice that 
the one moving right-handedly about the upward pointing axis is displaced above the center 
sitting hke a halo about it while that moving left handedly is displaced below the center hke 
an Elizabethan ruff below the head. One might have supposed that for j 3> Q* these two 
circular orbits would approach the central plane but although the cone becomes much flatter 
and more open the displacement between the direct and retrograde orbits actually increases. 
For circular orbits at distance a from S we have, for a Newtonian potential, = GMarn^ 
and the displacement is 

2j • r = a/^GMamlQ-^ + 1 ^ m^^Q.^J a/ {GM). 

We have been led to the case rrioNr^ — Q* — constant for reasons of mathematical simplicity 
but this case is more than a mathematical curiosity because: 

1. Of all the forces of Newton's N type [see Eq. (1.1)] only those of the form — v x 
rr~^(5*(^, (p) derive from a Lagrangian. For a monopole Q* is constant. 
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2. We may rewrite this force in the form 



NL = -g^v X r/r^ = — v x 



(1.17) 



where 




Qr/r^ ; Q^Q^c/rriQ. 



(1.18) 



We have introduced the velocity of hght c to make the analogy with magnetic forces even 
more obvious. Bg is clearly the field of a magnetic monopole of strength Q but since this sort 
of magnetism acts not on moving charges but rather on moving masses; it is a gravomagnetic 
field. Such fields are well known in general relativity see Landau & Lifshitz Theory of Fields 
(1966) §89 problem 1. They are position dependent Coriolis forces associated with what 
relativists less helpfully call the dragging of inertial frames. The field B^ as we have defined 
it has the same dimensions as g the acceleration due to gravity and Q/G has the dimensions 
of mass. In electricity, like charges repel while in gravity, like masses attract. It is the 
same with like magnetic monopoles, they repel while the gravomagnetic monopoles of like 
sign attract one another, hence the negative sign in Eq. (1-18) is best left there rather 
than combined into a new definition of the pole strength Q. We may find the Lagrangian 
corresponding to the force (1.17) by analogy with the electro dynamic case. There we add 
a term qw ■ A/c where q is the charge and A is the vector potential. For any poloidal axi- 
symmetric magnetic field one may choose A to be of the form where is the azimuth 
around the axis. We require 



from which one readily finds A — Q{1 -\- cos 9) gives the right B^,. Thus a Lagrangian for 
Eq. (1.1) is 



Although the dynamical system is spherically symmetrical the Lagrangian is not and can 
not be made so. The only spherically symmetrical vector fields are /(r)r. If A were of this 
form its curl would be zero and therefore could not be the field of a monopole. Of course 



Qv/r'^ = Bg = V X {AVcfy = x V0 



(1.19) 



C = Imov^ - moV{r) + Q^{1 + cos 9)v ■ V0. 



(1.20) 
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we can choose any axis we like and measure 9 and appropriately from it. The A field 
will then be quite different but it will give the same field by construction. Thus the 
difference between any two such A fields will have zero curl showing that A' = A + Vx i.e., 
a gauge transformation. The Lagrangian (1.20) is neither spherically symmetrical nor gauge 
invariant but it is a member of a whole class of equivalent Lagrangians with different axes 
which are related by gauge transformations. Whereas none of these is individually spherical 
the class of all of them is spherically symmetric. The moral is that it can be restrictive to 
impose symmetry on a single member of the class if the member is not gauge invariant. 

So far everything holds for any spherical potential V{r). We could for example choose it 
to be Henon's (1959) isochrone potential 2aVo/(a + -\/'"^ + o?) or its better known limits the 
simple harmonic oscillator a ^ r or the Newtonian potential a <^r. For all isochrones the 
orbits can be solved using only trigonometric functions (see e.g., Lynden-Bell 1963, Evans 
et al. 1990). Here we shall stick to the Newtonian potential V/rriQ — —GM/r. We have 
already shown that the motion lies on a cone whose semi-angle is given by cos~^(Q*/|j|); 
furthermore if we slit that cone along = and flatten it, the orbit will be exactly what 
it would have been in the absence of N i.e., a conic section. Of course when we slit and 



flatten the orbit's cone a gap appears whose angle is 7 = 27r 1 — L/yL^ + , see Fig. 2. 
An ellipse with focus at S and apocentre at 99 = would get back to apocentre at (/? = 27r 
but unfortunately the gap intervenes. On the cone we identify </? = not with (f — 27r but 
rather with = 27r — 7. Thus on the cone the ellipse will precess forwards by an angle 7 in 
each radial period. Fig. 3. This angle 7 is an angle like (p measured at S within the cone's 
surface. It is perhaps more natural to measure angles 77 around the axis of the cone; these 
angles are related through f] — (p/ cosx = L/(mor^cosx) so 77 = (/7secx = v\^\/L. 
In these terms the precession per radial period is 

Ar7 = 27r(|j|/L-1). (1.21) 

Newton in his proposition on revolving orbits showed that the addition of an inverse cube 
force led to an orbit of exactly the same shape but traced relative to axes that rotate at a rate 
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proportional to in the original orbit. It is natural to ask whether such an additional force 
can stop the precession around the cone of an orbit in the monopolar problem and so yield 
an orbit that closes on itself in fixed axes. Wonderfully a simple change in V{r) does this not 
just for one orbit but for all orbits at once. We thus obtain a new superintegrable system 
in which all bound orbits close. By analogy with Hamilton's derivation of his eccentricity 
vector (Hamilton 1847) we take the cross product of the equation of motion (1.1) with 
j = L + (5*f. On the right hand side two terms are zero and the remaining two are multiples 
of dr/dt cf. (1.12) so we find 

moj X dh/df = -(morV + Qlr-^)dv/dt. (1.22) 

This will integrate vectorially if the bracket is constant. Calling it GMm^ we find the 
potential must be of the form 

Wmo = -GMr-^ + -^r"". (1.23) 

2 uiq 

Evidently the required inverse cube repulsive force is proportional to the square of the 
monopole moment Q. Integrating (1.22) we have 

dr/dtxj = GMmo{r + e) (1.24) 

where e is the vector constant of integration. Dotting (1.24) with r we have 

4/r = (l + e-f) (1.25) 

where — L ■ j/{GMml) — const. Equation (1.25) is the equation of a conic section of 
eccentricity e which defines the direction to peri centre. But we have not yet proved that the 
orbit lies in a plane so (1.25) actually defines a prolate spheroid, paraboloid or hyperboloid. 
Nick Manton, by analogy with his work on monopoles in Euchdean Taub Space (Gibbons 
& Manton 1986), showed us that the motion is in fact planar; for using (1.15) Eq. (1.25) 
becomes on multiplication by Q^r 

QJ* = (j + <5*e) • r 
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which demonstrates that the orbit hes on a plane whose normal is j + Q^e. As r also lies on 
a cone this provides another proof that the motion lies along a conic section. 

Notice that the vector integral e in (1.24) together with the integral j appears to provide 
six integrals of the motion. However they are not all independent because — e-j = f • j = Q*. 
So they provide 5 independent integrals. Thus we have a new superintegrable dynamical 
system in which the bound orbits exactly close (cf. Evans 1990, 91). 

It was the beauty and simplicity of these results for monopoles in classical mechanics 
that led us to believe that a similar simplicity might well be discernible both in quantum 
mechanics and in general relativity. We were not disappointed, both had already attracted 
attention. Hautot 1972 discusses the separation of variables in r, 9, coordinates. The 
vector integral j is preferable because the generahty of motion on cones is then seen. For 
motion in special relativity j is still conserved provided L is interpreted as mor x dr/dr = 
mor X v/ ^JT^^vy^. Goddard and Ohve (1978) in their excellent review of monopoles in 
gauge field theories quote Poincare (1895) for this integral in the classical case of a pure 
electromagnetic monopole. 

II. DIRAC'S MONOPOLE AND THE SPECTRA OF MONOPOLAR ATOMS 
A. Gauge transformations, Schrodinger's equation &; Dirac's quantised monopole 

The Lagrangian for a particle of mass mo and charge — e in an electromagnetic field is 

>C = |mor^ — er • A/c + e$(r) 

where $ is the electrostatic potential and B = Curl A is the magnetic field. The momentum 
conjugate to r is 

p = dC/df — mor — eA/c 

which is not a gauge invariant quantity. However the particle's momentum moi" = p + eA/c 
is gauge invariant and therefore has greater physical significance. The Hamiltonian is given 
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by 



= p • V - £ = {2mo)-\p + eA/cf - e$. 



(2.1) 



For Schrodinger's equation we replace p by —ihV and solve Hijj — Eijj for the wave function 
of a steady state. A given magnetic field B can be described by many different vector 
potentials A related by gauge transformations A' = A + Vx. Each will give us a different 
Hamiltonian. Let us first see how the different wave functions corresponding to these are 
related. Define a new function ip' such that 



then {—ihV + eA/c)tjj — exp[iex/ {hc)]{—ihV + eA/c+ eVx/c)i'' and the combination 
A' = A + Vx has appeared. Applying the above operator twice we see that 



where H' is H with A replaced by A'. It follows that Schrodinger's equation Hip — Eip 
implies H'lp' = Eip' so under gauge transformation ip transforms to ip' given by (2.2). 

This is a perfectly good wave function whenever x is single valued but following Aharonov 
and Bohm (1959) we now consider the wave function of a particle outside a small impen- 
etrable cylinder R = a. If we take Vx = V{F(j)/27r), R > a where ([) is the azimuth this 
corresponds to the same magnetic field outside the cylinder but a different magnetic flux 
within it because 



which identifies the constant F as the extra flux threading the cylinder. If we adopt our 
transformation of wave function for a gauge transformation we get the phase factor 



ip — exp[iex/(^c)]'0' 



(2.2) 



Hi/j = exp[zex/(^c)]//V' 




= B-dS + F 



exp[— ieF0/(/ic)] 



(2.3) 



which is only single valued when the flux takes the special values 



F = N{hc/e) 



(2.4) 
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where N is an integer (positive, negative or zero). (This N is not the force coefficient of 
Section I). Thus while we get the correct wave function for those particular values of F, we 
need to solve the problem anew with the correct boundary condition that ip' must be periodic 
in whenever F is not an integer multiple of the flux quantum hc/e. Indeed when it is not, 
there is interference between the two parts of a beam of electrons that pass on either side of 
such a cylinder just because their phases differ by e § AA ■ dl/{hc) = eF/{hc). It was just 
this phase shift that was observed in the experiments demonstrating the Aharanov Bohm 
effect of the magnetic flux even when the electron beams were untouched by the magnetic 
field. There is an intimate connection of this result with Dirac's (1931) earlier quantum of 
magnetic monopole from which one flux unit (2.4) emanates.. This comes about because in 
the presence of a monopole ^ A ■ d£ is itself multivalued. 

Consider the integral § A - dl around a small loop; this is clearly the flux of B through 
the loop but such a flux is ambiguous in the presence of a monopole since it depends on 
whether the surface spanning the loop is chosen to pass above or below the monopole i.e.. 
Si or 5*2 in Fig. 4. The difference between these two estimates is just js^_s^ B • = AnQ 
by Gauss's theorem. Inserting this A / A • c?£ in place of F in (2.3) we see that the wave 
function will only have an unambiguous phase provided 

47rQ = N{hc/e) (2.5) 

i.e., provided that the monopole strength is quantized in Dirac units of \Tic/e ~ 

The quantum of magnetic flux (2.4) is inversely proportional to the charge. Quanta of 

half this size are observed in the Josephson effect in superconductivity where the effect is due 

to paired electrons. There is some evidence for the larger unit (2.4) in ordinary conductors 

at low temperatures (Umbach et al. 1986). 

Returning to Schrodinger's equation (2.1) and using the vector potential [cf. under 

(1.19)] 

A = -Q(l + cos^)V(/) (2.6) 
we have the correct magnetic fleld for a monopole of strength Q but we notice that A is 
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singular along the line ^ = although it is regular along 9 — n. Near the singular hne 
A — 2(5V0 which is the vector potential of a tube carrying a flux AttQ downwards, 
thus formula (2.6) represents the vector potential of a magnetic monopole fed its flux by the 
singular half hne ^ = 0. This half hne gives an unobservable Aharanov-Bohm effect provided 
AnQ — Nhc/e that is provided the monopole is a multiple of the Dirac (1931, 48) unit. Extra 
interest in his monopole comes from his argument that it can also be read as a reason for 
charge quantization, because, if Q is the least monopole, then e must be a multiple of he/ Q; 
thus in his picture, charge quantization and monopole quantization spring from the same 
source. It is of interest that A in (2.6) is single valued. It has avoided the multi-valuedness 
alluded to above by having the singular string at ^ = down to the monopole. This plays 
the role of the cut in multivalued functions in the complex plane Wat & Yang (1976). An 
interesting historical remark is that Schrodinger in 1922 saw that quantum conditions in the 
old quantum theory led to F = ^[f ^dt — A • dx] — nh while Weyl's gauge theory led him to 
consider exp (— r/7) with 7 as yet unspecifled. He realized that the identiflcation 7 = —ih 
would lead naturally to such quantum numbers and after de Broglie (1925), he built on this 
idea to invent his wave mechanics in 1926. (See Yang 1987). 



B. Solution of Schrodinger's Equation 



Written in spherical polar coordinates Schrodinger's equation is 



d 



+ 



2mor^ [ dr 



d 

djjL 



'1 2^^^ 



+ 



Q^2 ' -'Q^ 4 



(2.7) 



here ^ has been written for cos 6 and N is the number of Dirac monopoles on the nucleus. 
(p only occurs as d/d(j) in the above equation so we may take one Fourier component with 
■0 oc exp [imcj)) and m an integer positive, negative or zero in order that if) be single valued. 
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On multiplication by — 2mor^/ {h'^ip) we then find the separated equation 



ip d/j, 



[m-iVi(/x + l)p 



2 



l|_/,^\_r^ (2.8) 



Writing %l) — ipr{r)ipn{ii) the left hand side is a function of /i alone and the right hand side is 
a function of r alone so both must be a constant which we call — C. The resultant equation 
for ipf_i has regular singular points at ^ = ±1. The indicial equations for the series solutions 
about fjL — ±1 have regular solutions behaving as (1 — //js'"*"^! and (l + /x)2l"*' respectively, 
so we remove those factors by writing 

V'^=(l-/.)il"^-^l(l + /.)iHi.(^). (2.9) 

After some algebra the equation for F takes the form 

(1 - ii^)F" + [(|m| + - ii) - {\m - N\ + + ii)]F'+ 
+l[2C - m{m - N)- \m\\m - - |m - iV| - |m|]F = 0. 

(2.10) 

We now write z — |(1 + jj), so z{l — z) — {1 — /i^)/4 and dz — ^dji which reduces the 
above equation into the standard form for the hypergeometric equation i.e., 

z{l - z)d^F/dz^ + [c - (a + 6 + l)z\ dF/dz - abF ^ (2.11) 

where 

c=\m\ + l (2.12) 

a + b^\m\ + \m- N\ + l (2.13) 

and —2ab is the final square bracket in Eq. (2.10). 

The hypergeometric function finite at /j, — —1, z — diverges hke (1 — jjiY~"'~^ at = 1, 
that is twice as fast as the first factor in (2.9) converges, so in order to get convergence the 
hypergeometric series must terminate. This occurs only if a or 6 is a negative integer or 
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zero, w.l.g. taking b — —K we find that F reduces to a Jacobi polynomial P^^(//) so that 
i/jfj, takes the form 

V'M = C,mn{^ - + ^^Mn^l p\^-NUm\ ^2.14) 

here J^^ '^fid/J- — 1 and Ckmn is the normalization 

1 



{2K + 


m- N\ + 


m 


+ 1)K\{K + 


m — 


N\ + 


m\)\ 


2|m-Ar|+|m|+l 


{K + 


m-N\)l{K + 


m\y. 



The condition that b — —K gives 

-2ab = 2K{\m\ + \m-N\ + l + K) 

and hence (noticing that K — leaves (2.14) finite) we have 

C = K(K + 1) + K(\m\ + \m- N\)+ 
+ ^[m{m — N) + |m||m — A^l + |m — + \m\] . (2-15) 

If we write j = K + |(|m| + \m — N\), then we notice that j is a positive half-integer and 
j > l{\m\ + \m- N\) 

C^j{j + l)-NyA. (2.16) 

Thus C and j are only integers when N is an even integer. When is odd C and j are an 
integer ±|. For given j and N > 0, rn — ^ takes the 2j + 1 values from —j to +j in steps 
of 1. For N — 1 the ground state has j = \ and C — \ rather than the values familiar 
from the normal hydrogen atom. The j — \ states with m — 1 and m = are degenerate, 
see Fig. 5. 

With the value (2.16) for C we now turn to the radial equation for ipj., (2.8). Here the 
treatment is very close to the classical case clearly laid out by Pauling and Wilson (1935). 
We take ^ — Ze/r, Ze being the nuclear charge, and E negative. We write 

= -2mQElh^ (2.17) 

C = m^Ze^h-'^ar^ (2.18) 
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and use a normalized radius f = 2ar. As all the radii in the rest of this section are so 
normalized we shall forget the ~ and take it as read. Eq. (2.8) now takes the form 

The asymptotic form of this equation for large r is ■0r = ■0r-/4 so ■0r- ^ exp(r/2) or exp — 
(r/2). Of these only the second is acceptable so we write 

i^r = exp(-r/2)r7(r) (2.20) 

where /(r) may be expanded in series about the origin in the form 

oo 
p=0 

and s is chosen so that Oq 7^ 0. The indicial equation found by substitution of the series 
(2.20) into (2.19) is 

[s{s + 1) - C]ao = 
but by hypothesis ao 7^ so using (2.16) s is given by 

(5 + I)' = i + « + 1) = I + C = (i + 0' - \N\ (2.21) 

The recurrence relation for general p is then 

p{p + 2s + l)ap = (s + p - C)ap-i 

and the asymptotic form for large p is ^ ap-i/p which shows that f ^ e^. In that case 
ipr would diverge at large r. This is unacceptable so the series must terminate. Thus there 
must be a positive integer p = n' + 1 such that 

( = p + s^n' + s + l (2.22) 

with s given by (2.21). Returning to (2.18) and (2.17) this gives the eigenvalues for the 
energy in the form 

p_ rnpZ^e^ 1 _ mpZ^e' 1 

2h^ {n' + s + l)^ 2h^ {n + Ay ^ ' ' 
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where n — n' + J + 1 where J takes values 0, 1, 2 ... replaces the usual i and 

J = >o 

N is the number of Dirac monopoles on the nucleus and 

A = V(-^ + i)(J+| + |A^|)-(j+|) >0. 
Notice that A depends on J as well as jA?"! and is only zero when N — 0. For large J/\N 

A^l\N\ 

while for the ground state J — 

A = i(^2|7V| + i-i) 

which becomes | (^/3 — 1^ ior N — 1. So A is not small. For a spinless electron the 
degeneracy of a state of given J and nis2j + l = 2J + l + |A^| with m — ^\N\ taking values 
from —j to +j. Notice that the ground state J = 0, n = 1 is a doublet for = 1 and has j 
value I with m — and m — +1 states even before we have allowed for further degeneracy 
due to electron spin. A single Dirac monopole thus gives some effects reminiscent of spin | 
particles (Goldhaber 1976). 

The dependence of A upon J lifts the degeneracy of the different J states {i states) that 
occurs in normal hydrogen. The energy levels are near to those for an atom with a true 
spinning electron laid out in Figs. 6, 7, 8 and Tables 1 & 2. The degeneracy would return if 
the extra repulsive potential |(5^/(mor^c^) were included. Then the —\N'^ in (2.21) would 
be cancelled so s would become equal to j. 

C. Selection Rules 

The string to the monopole makes it look non-spherical but this is not truly the case as 
putting the string in any other direction can be achieved by a mere gauge transformation. 



^ 1 \N\ 1 / |Ar| y 

4j+i 8[j+y 
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Therefore without loss of generality we may evaluate transition moments by taking the 
displacement in the z direction in which case we get 

The radial integral is that for normal hydrogen but the scales have changed since s in (2.22) 
is no longer I but is given instead by (2.21). We shall concentrate on the important change 
in selection rules given by the final integral. 

Whereas for the Legendre Polynomials in normal hydrogen wave functions we have 

so that f^i Pe'uPedfj, is only non-zero when £' — £ — ±1, for the Jacobi polynomials in 
monopolar hydrogen 

l,Pf = (a^Pf^, + a^Pf + a,Pf_,) /a^ 

where 

ai^2(K+l)(K + a + (3 + l){2K + a + (3) 
a2 = {2K + a + p + l){a^ - p^) 

{2K + a + P){2K + a + P + 1){2K + a + P + 2) 
04 = 2{K + a){K + (3){2K + a + P + 2) 

so that 1+^(1 - //)"(1 + nfP^^iiPK^dii will be non-zero when K' - K ^ ±1 or 0. [The 
term is only absent when a2 = i.e., a= \m — N\ = j3 = \m\. This occurs for iV = always, 
for N — 2 when m — 1, but never for N — 1.] 

Thus there is a significant change in the selection rules for electric dipole transitions. 
Some might imagine that magnetic dipole transitions should be important but the mag- 
netic monopole is on a heavy nucleus and barely responds to an oscillating magnetic field 
so it is still the electric dipole transitions of the electron that are important, m is un- 
changed for a dipole along the z-axis so AK = ±1 or leads directly to Aj and hence 
AJ = ±1 or for such transitions. 
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Even order of magnitude estimates show that the interaction of the electron spin's mag- 
netic moment with unit monopole gives not a dehcate fine structure but significant changes 
in the eigenvalues! Thus to find the true eigenvalues the Dirac equation is a necessity! Be- 
fore treating it we clear up some details. We took the form (2.2) for A corresponding to 
a monopole with a string along /j, — +1. For |A^| > 2 we could have taken two or more 
inwardly directed strings of flux. Are these different string configurations really different 
monopoles or do they all give the same eigenvalues? The effect of such a change is to add a 
unit flux string along the z axis. It is simple to show that this is equivalent to adding one 
to m everywhere that it occurs. Provided we do that also to m in the definition of j under 
(2.15) the final spectrum remains unchanged. What does change are the K and m values 
associated with a given j value. 

A second detail is the value of mo which for A = would be the reduced mass of the 
electron so for hydrogen it is mo = memp/(me -l- m^). 

Particle physicists expect a heavy mass for any monopole so any monopolar hydrogen 
will have a nucleus much heavier than the proton and rUe should be substituted for mo in 
predicting spectra. A third detail for later reference is the energy spectrum of the relativistic 
Klein-Gordon equation. Here we follow Schiff 's treatment and obtain writing — Ze^/ (he) 



Returning to the classical conserved quantity cf. (1.14) we see the conserved quantity is 
not the particle's angular momentum L = r x mov but rather that supplemented by eQc~^r. 
The physics behind this supplement lies in the Poynting vector of the electromagnetic field 
which carries an angular momentum 




(2.25) 



where 




(2.26) 



D. Angular Momentum 
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}- fr'x (Ex^r')d\'^^ f (E ■ V) rW = 
4:TTcJ V r'2 / AttcJ 

~^ J r'V • Ed'r' = = +i7V;if 



47rc 

where V • E = — e47r(5^(r' — r). The total angular momentum is thus j = L + eQc~^r. 

As we saw above (2.1) mov = p + eA/c in the presence of a magnetic field 
so the operator representing j is r x {—ihV + eA/c) + jNhr. The commutators 
[—ihdj + eAj/c,—ihdk + eAk/c] — —ihec'^SjkiBi — —iheQc~^ejkix^ /r^ and [—ihdj + 
eAj/c,x^] — —ihdj enable one to derive the commutator 

[jjjk] = i^^jkdi 

which demonstrates that j obeys the angular-momentum algebra of the rotation group. One 
may also demonstrate that commutes with j and that j± = ± ijy are the raising 
and the lowering operators for j^. Prom which it follows by the usual argument that the 
eigenvalues of jz are —jh to +jh and that the eigenvalues of p are j{j + l)h^. But |jp = 
|Lp+iiV^/i^ so the eigenvalues of |Lp are j {j + 1) - jN'^l h^. Now looking at our separation 
of variables expression (2.8) we see that the LHS is just — ^~^|L^| by construction and hence 
C — + 1) — |-^^) which agrees with (2.16) and identifies the j defined there with the 
generalized angular momentum eigenvalue defined in this section. Note that for a single 
Dirac monopole and a non-spinning electron we showed (2.16) that j took half odd integer 
values. 

In the next section we look at the Dirac equation for a spinning electron. There the 
correct generalization is j = L -|- ^Nhr + ^ha. 

This new j obeys the angular momentum algebra of the rotation group but now its 
eigenvalues are j{j + l)ti^ with j taking integer (or half integer) values > depending on 
whether N is odd or even. 
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E. Dirac Equation 



The Dirac equation may be written in standard notation 



Hip = -ca • (p + eA/c) - pnioc^ + V ip ^ Eip 



with the newly defined j we find ^ = [j, i7] = so that each component of this generahzed 
j commutes with the Hamiltonian always provided that A is a vector potential for the 
monopole. Following Schiff's treatment (1955) we define Pr — r~^(r ■ p — ih) and ar — 
r~^{oi ■ r) and Kk = Pa' ■ (r x (p + ^) +hj- No A term is needed in pr since r • A = for 
our monopole. 

The Hamiltonian is now written 



and as before ar, (3 and Pr all commute with k. The eigenvalues of k follow by squaring the 
definition 



for this change of k the usual separation of the Dirac equation goes through unscathed and 
following Schiff one obtains the energy levels 



H — —carPr — i — o-rPk — PmoC^ + V 
r 





where s — {k"^ — a^) and = Ze^ / (he) . 

n' is the radial quantum number. Inserting our eigenvalues k'^ — (j + | 




kN^ with 



j ^ J+ i(|iV| + 1) and J = 0, 1, 2, etc. we have 




where n — n' + J + 1 and 



A = ^(J+l)(J+l + |iV|) 



-al-{J+l). 
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These energy levels were first derived by Hautot (1972, 73) generalizing Harish-Chandra's 
(1968) separation of the Dirac equation^. For the scattering by monopoles see also Gold- 
haber (1965), Kazama & Yang (1976) and Kazama et al. (1977). We have drawn the 
bound energy levels that result Figs. 6, 7 and 8 and derived the wavelengths of the lines of 
"Monopolar Hydrogen" with one or two Dirac monopoles attached to the nucleus. Tables 1 
& 2. Schwinger (1966) has suggested that the unit monopole should have the strength of two 
Dirac monopoles. With colleagues he has also calculated the motions of charged monopoles, 
dyons, under their mutual attraction (Schwinger et al. 1976). While monopoles may seem 
esoteric it is worthwhile looking for lines of monopolar hydrogen in the spectra of exotic 
astronomical objects. 

III. GRAVOMAGNETIC MONOPOLES IN GENERAL RELATIVITY, NUT 

SPACE 

A. NUT space the general spherically symmetric gravity field 

Zelmanov (1956) and Landau & Lifshitz (1966) in developing their very physical approach 
to general relativity consider stationary space-times and put the metric in the form 

ds'^ = e-2^(rfx° - A^dx'^Y - -f^^pdx'^dx^ (3.1) 

where i' > 0,Aa and 7^/3 are independent of = ct. (Our z/ is — |z/ of Landau & Lifshitz). 

However this form is not unique since a transformation of time zero x'^ = x^ + x(^") 
leads to 

ds^ = e-^^idx"^ - A'Jx'^f - j^pdx^'dx^ 

where A'^ = Aa + VaX so under such a change A undergoes a gauge transformation. Landau 
& Lifshitz also show that 7^^ can be regarded as a metric of space i.e., the quotient space 

^For the further generalization to the problem with an additional Bohm-Aharonov string, see 
Villalba (1994, 95) & Hoaiig et al. (1992). 
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V^/V- (Geroch 1970) - as opposed to space-time. They show that test bodies following 
geodesies of space-time depart from the geodesies of space as if acted on by gravitational 
forces which in our notation take the form 



f = 



mo 



E„ + - X e-^B„ 
. c 



where the gravitational field 



and 



Eg = Curl A. 



(3.2) 



(3.3) 



(3.4) 



The conserved energy of the particle in motion is 



e = moc^e"'' (l - v'^/c^) 



-1/2 



(3.5) 



e < 1 is the redshift factor by which energy is degraded. 

Rewriting Landau & Lifshitz's form of Einstein's equations (§95 problem) we find 



div Bg = 



Curl Eg = 



(3.6) 



(3.7) 



div Eo = -c 



-2 



47rG 



1 — v"^/ (? 



2^ - E^ 



(3.8) 



where p is the energy density in the rest frame of the fluid, 3p is the trace of its pressure 
tensor and v its velocity defined locally by local time synchronized along the fluid's motion. 
For non-relativistic velocities this equation reduces to Poisson's equation with the primary 
term on the right being AnGp. The remaining term has the form of a negative energy density 
contributed by the gravity fields. The next equation takes the form 



Curl (e-^'Bg) = -c-^[16nGjg - 2cEg x e-'^Bg 



(3.9) 
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Notice that e'^'Bg occurs also in that combination in the expression for the force. It is 
attractive to regard the final term as an energy current corresponding to a Poynting vector 
flux of gravitational field energy, jg the matter energy current is given by 

. _ pc^ +p 

The final Landau & Lifshitz equation for the 3 stress tensor is 

P"/3 _ — _|_ ij^^oiP^ _|_ 
+ e-^^iB^B^ - 5^7"^) + E^E^, (3.10) 

P°'^ is the 3 dimensional Ricci Tensor constructed from the metric 7"^. Those familiar with 
the Maxwell stresses of magnetic and electric fields in say magnetohydrodynamics will find 
some interest in the field terms on the right. The matter terms may be rewritten as physical 
quantities for an isotropic fluid in motion 

It should be stressed that all these equations hold good even when space-time is strongly 
curved. Unlike some treatments they are not restricted to nearly flat-space but it is assumed 
that the space-time is stationary. 

To flnd the general spherically symmetric solution for empty space we take dl'^ — 
'^apdx°'dx^ — e'^^dr'^+r'^dP where f is the unit Cartesian vector (sin 9 cos 0, sin 9 sin 0, cos 9). 
Then dr'^ — d9^ + sin^ 9 dcj)^ but the advantage of the vector notation is that no axis for 9,(f) 
need be taken. In spherical symmetry Bg must be radial and divergenceless so Gauss's 
theorem gives |Bg|r^ = Q — const which is the fleld of a gravomagnetic monopole 

Bl = -Qe-'/r'. (3.11) 

Reinserting = c^Vu into (3.8) we have 

= -u" + u'^ - 2u'/r + AV + ^e^^^-''^ {cr)'^ = 0. (3.12) 
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rpaP _|_ Irp^a/S _ 



SttG 



To form we need the 3 dimensional Christofel symbols 



= + 7,.,M - 7m.,.) (3.13) 

which are 



(3.14) 

(3.6), (3.7) and (3.9) are identically satisfied. The surviving equations of (3.10) are 

Rrr ^ + j,'2 ^ _ 3^/^ ^ Q (3.15) 

and 

r r 2 

+ i/'e-^^ = 0. (3.16) 

Equations (3.12), (3.15) and (3.16) must be solved for u and A. Eliminating u" from (3.12) 
and (3.15) we find 

2(A' - v') + \e^^^-'^Q^c-\--'' = (3.17) 
which integrates on division by e'^^^"'^^ giving 

e-2(A-.)^_gV2 + C (3.18) 

where q — Q/2c^ which has the dimensions of a length. Multiplying (3.16) by e^^ and using 
(3.17) and (3.18) we have 

(C - gV-2)(r-i - 2u') = gV'^ + e+^^/r 
dividing by e^'^(C — q^r'"^) we obtain 

(e ) + - ^r;^ — IT e - — = 
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which is hnear in e and readily solved by integrating factor to give 

= ^ - 1$ + ^V^;^ (3-19) 
where C and C are both constants. It follows from (3.18) 

7,, = e^^ = (C - gV2)-ie+^^ (3.20) 

To get e"^'' and 7q/3 asymptotically of Schwarzschild form we need C — 1 and C = — m, the 
asymptotic mass GM/c^. Thus we find 

= e-''^ = 1 - 2r-' {(f + mVr^ - P ) (3.21) 

7.. = (1 - ?V-2)e+^'^ (3.22) 

which are the metric components of NUT space. Notice that when Q — 2qc^ — this reduces 
to Schwarzschild's metric. The metric is completed by taking a vector potential Aa for the 
gravomagnetic field B^; any one will do since they are connected by gauge transformation 
which merely changes the zero point of time. As we saw in Section I it is impossible to choose 
a spherically symmetric vector potential but this does not affect the spherical symmetry of 
the physics. A suitable vector potential is that given in (1.19) which gives us the metric 

ds^ = e-^" {cdt - 2q{l + cos9)d(l)f - 
-(l-?Vr2)e+V-rW (3.23) 

where e~^'' is given by (3.21). This metric is more commonly written in terms of the radial 
variable f = \^r^ — q'^ because the square roots disappear leaving an analytic expression, 
however we have preferred the variable that makes the surface area of the sphere 47rr^ as in 
Schwarzschild space. Of course the metric (3.23) appears to have a preferred axis but this is 
illusory because we can switch it into any direction we like by a gauge transformation, see the 
discussion under (1.20). The horizon where goo changes sign is given hy f — m + \/q^ + rn? 
at which point ^rr changes sign also as in Schwarzschild space. 
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NUT-Space was discovered in Ehlers' thesis (1957) and rediscovered by Newman, Tam- 
burino & Unti (1963). It is closely related to Taub's (1951) metric and their relationship 
has been beautifully illuminated by Misner & Taub (1969). The fact that NUT space has 
a gravomagnetic monopole was found by Demianski & Newman (1966), who also found a 
NUT version of Kerr space. See also Dowker & Roche (1967). 



B. Orbits and Gravitational Lensing by NUT Space 

The geodesies of NUT space may be determined from S J ds — using the metric in the 
form (3.1). When ds'^ 7^ we write r for the proper time and when ds'^ 7^ we replace it by 
an affine parameter (which we also call r). Varying t = dt/dr and using the fact that the 
metric is stationary we have 



e {ct — Aax"") — e — constant. 



(3.24) 



Varying we find 



+ 



1 d 



e-^''{ci-Apx^)A„ + j„px^ 



+ 



(3.25) 



Using (3.24) and transferring the edAa/dr — ex^^dpAa term into the second bracket we find 
the equation of motion in which A only occurs through daAp — dpAa — rja/SyB'^ where rjapj 
is the antisymmetric tensor, times the alternating symbol. 



dr 



21 Pi ,oc^^ ^Va/S-f 



dx' 
B^x^ 



= 



(3.26) 



We now write 7q,/j in the form involving the unit Cartesian vector r. 



^ocpdx'^dx'^ = e^^dr'^ + r'^{drf. 
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5r the variation of f is an arbitrary small vector perpendicular to f . Thus making variations 
with r fixed we deduce from (3.26) Using (3.11) for 







since 5v is an arbitrary vector perpendicular to f we deduce that 



fx— r — 



dt \ dr 



) 



dL 
d^ 



d 

d^ 



{eQr) 



L + eQr — j— const. 



(3.27) 



Except for the factor e which reduces to moc^ in the non-relativistic case we see that this is 
precisely the vector integral (1.14). Dotting Eq. (3.27) with f we find j ■ r = eQ showing 
that f lies on a cone similarly L • j = = — s^Q^ = const so L moves around a cone. 
The radial equation of motion is redundant since we may use the energy and the equation 
[ds/drY — U — 1 or instead. C/ is 1 for time hke geodesies and for hght-like ones. 
This gives us 



To see the geometry of the trajectory we introduce the curvilinear angle if of §1 measured 
around the cone's surface. Then r^<^ = L so Eq. (3.1) can be integrated by quadrature 



In general this integral can not be performed explicitly for the A and i/ of NUT space 
even after substitution in terms of f to make it more analytic. We therefore turn to the 
^ + rn? limit well away from the event horizon. This is the important case in all 
gravitational lenses observed to date. In that hmit the /r^ term in the effective potential 
is attractive and therefore of the wrong sign to give the non-precessing orbits of Section I. 
The precession around the cones is faster than in the classical Kepler + monopole problem 
by a factor 3/2. To the first order in m/6 where h is the impact parameter, we find a 
bending angle measured like </? of Ac/? = 4m/6 just as for the Schwarzschild metric; however 




r2 -2 

— L r — 



U. 



(3.28) 




(3.29) 
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the difference is that this angle is measured around a cone not in a plane. Again to first order 
we can find the effect of the gravomagnetic field by integrating the momentum transfer along 
the unperturbed straight line path. This gives an out-of-plane bending of 4g/6; a result that 
is confirmed by the full NUT space calculation. Nouri-Zonoz, M. & Lynden-Bell, D. (1997). 
Thus the major effect of the gravomagnetic monopole Q is to twist the rays that pass it. 
While the bending angle is proportional to the effect is exaggerated when looking down 
the line toward the NUT lens by the factor Di/h, so the twist around the lens is AqDi/b^. 
Here D^is the distance from the observer to the lens. The same exaggeration factor occurs 
for the normal gravitational bending so for a source at infinity and an image at (6, 9) in the 
plane of the sky at the lens's distance, the apparent position of the source is 



This expression defines a map from image to source. Prom this map one can work out both 
the shear and the magnification of a NUT lens in the large impact parameter regime. The 
magnification of area and thus luminosity is 

db''/dbl = [1 - m-^Dl{m^ + q^)]-\ (3.30) 

A small circular source will be imaged into an ellipse of axial ratio 

6^ + 4 DL{m + VmFTq^) 
62 + 4 Di(m - Vm^TF) 

with the short axis of the ellipse inclined to the radius at the angle (see Fig. 9) 

tan~^ 



+ \Jm^ + q^ ) 

This is 45° for q ^ m and 13° for Q = 2qc^ = mc^. This spiral conformation of the 
images about a NUT lens is very characteristic. It is not displayed in normal gravitational 
imaging and the gravomagnetic lens due to a rotating object seen pole on does not show it 
because the twist of the ray as it approaches such a lens is cancelled by the opposite twist 
as it recedes. Thus the discovery of a spiral shear field about a lens would indicate the 
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presence of a gravomagnetic monopole. Such effects should be looked for by those studying 
gravitational lenses. The expectation must be small but the reward might be an amazing 
discovery. 

C. Quantization of Gravomagnetic Monopoles and their Classical Physics 

By analogy with Dirac's argument for the quantization of magnetic monopoles and 
charges, Dowker & Roche (1967), Dowker (1974), Hawking (1979) and Zee (1985) have 
suggested quantization of gravomagnetic monopoles and energy. Corresponding to Dirac's 
Qm^ = \NfLC for magnetic Qm, they have Qruo = \Nhc for gravomagnetic monopole Q. 
This implies that both Q and mass mo are quantized in conjugate units Qi and mi obeying 
Qirrii = ^hc. Whereas such ideas are naturally attractive they do not naturally lead to 
a self-consistent relativistic theory. For instance, looking at the Klein-Gordon equation in 
NUT space and separating variables with ^jj oc e*('"'^+'^*), one finds an eigenvalue equation 
for uj. The Dirac monopole quantization condition, Q{Tiuj/c^) = ^Nhc with an integer, 
shows us that the only possible eigenvalues u are integer multiples of ^c^/Q and the cor- 
responding energy hu is the total energy of the 'orbit' including rest mass. However this 
condition conflicts with the energies of the bound states^ which are not integer multiples 
of any unit. Mueller & Perry (1986). Thus if such ideas are viable at all a more radical 
change in basic theory is needed. In + + -{ — NUT space it does not appear to be possible to 
build a consistent quantum theory like Dirac's magnetic monopole theory. This is what Ross 
(1983) concluded and is related to Misner's (1963) finding that NUT space contains closed 
timelike lines, with time being periodic every Sirq/c. For a discussion of the energy levels 
in + + ++ Taub-NUT space, the reader is referred to the papers by Gibbons & Manton 
(1986). This space was shown to be relevant to the interactions of monopoles by Atiyah 

^To get definite bound states one must impose a potential barrier so that the black hole is not 
reached. 
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& Hitchin (1985). §If magnetic monopoles exist, Maxwell's equations must be changed to 
include div B = 47rpm, Curl E + = 47rjm, where pm is the monopole density and jm is 
the monopole current density. Such modified Maxwell equations do not come with a vector 
potential. It is natural to ask how general relativity must be modified to allow for gravo- 
magnetic monopole densities and currents. While this is not so obvious we conjecture the 
generalization will be to spaces with unsymmetric affine connections which have non-zero 
torsion. It would be interesting to demonstrate this conjecture as it could introduce a greater 
degree of physical understanding of those spaces. 

IV. OBSERVABILITY 

Following Kibble's (1980) suggestion that magnetic monopoles would be a natural con- 
sequence of the Big Bang, they have long been sought. 

We have concentrated on the spectra of monopolar atoms and the lensing properties of 
gravomagnetic monopoles since these are ways in which, at least in principle, monopoles 
might be discovered observationally. Spectroscopically one may argue that the best place 
to look is in the spectra of supernovae, quasars or active galactic nuclei where the basic Ly 
a lines of Table I or II might be seen as very weak absorption lines in very high resolution 
spectra. Quasars have the advantage that these lines will be shifted into the visible. We have 
looked at lUE ultra-violet spectra of Supernova 1987A and seen no lines at the wavelengths 
2774.62 or 2733.78. More supernovae and stacked high resolution spectra of quasars should 
be pursued. Although in regions of observed magnetic fields the limits obtained spectro- 
scopically will fall far short of the Parker (1970) bound. While the nature of the dark matter 
that constitutes most mass in the universe remains unknown, such esoteric possibilities are 
worth pursuit. 

Searches on Earth have produced one unrepeatable event and a monopolar observatory 
under the Grand Sasso that has so far found no monopoles in Cosmic Rays. There has been 
a speculative suggestion, Kephart & Weiler (1996), that the leveling up on the numbers of 
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cosmic rays at the highest energies might be due to monopoles but there is no confirmation 
of that idea. To date the best hmit on the numbers of monopoles in interstellar space 
comes from the Parker (1970-71) bound. This arises from the idea that too many magnetic 
monopoles would 'short out' the galactic magnetic fields that are observed. A good general 
discussion of such limits may be found in the book of Kolb & Turner (1991). For more 
recent work on monopoles in field theory see reviews by Olive (1996, 97), and the papers 
by Sen (1994) and by Seiberg & Witten (1994). More details of the fundamental work on 
monopoles in field theory by 't Hooft (1974) and by Polyakov (1974) can be found in the 
review by Goddard & Olive. 
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FIGURES 

FIG. 1. The circular orbits about a central potential endowed with a monopole. The orbits in 

opposite senses are displaced above and below the centre of force. For a Newtonian potential their 
vertical separation gradually increases as their radii are increased. The orbits with a given angular 
momentum |L| lie on cones with opening angle cos~-'^(Q*/|j|). 

FIG. 2. When one of the cones is slit and flattened a gap opens along the slit. On the cone 
itself the sides of this gap are identified. Orbits which close on a plane will not close on the cone 
because of the gap. As a result they precess. 

FIG. 3. An ellipse precessing around a cone of semi angle 70° making a rosette orbit on it. 

FIG. 4. A monopole and its B field showing the surfaces 81,82 and 83 = 81 — 82- 

FIG. 5. j values allowed by the conditions j > ^{\m\ + |m — ^|). j cannot be less than the 
average of the two faint V lines in the diagram. 

FIG. 6. Energy levels for a spinning electron in hydrogen with 0, 1, 2 or 3 Dirac monopoles on 
its nucleus. Excepting 'isotopic' shifts due to the changed nuclear mass and relativistic corrections, 
the energy levels of the ground states are in the ratio 1 : ^ : | : ^. 

FIG. 7. Energy level diagram E{n, J) for N = 1, hydrogen with one Dirac monopole on its 
nucleus. The nucleus has been assumed to be fixed. 

FIG. 8. Energy level diagram E{n, J) for N = 2, hydrogen with two Dirac monopoles on its 
nucleus. 

FIG. 9. Gravitational Lensing by NUT space of a small circular source at 8 appears as an 
inclined ellipse at the image /. Many such images make a spiral effect around the NUT lens L. 
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TABLES 

TABLE I. Wavelengths in Angstroms of the two Lyman Series, the five Balmer Series and the 

eight Paschcn Scries of hydrogen with (N = 1) one Dirac monopole attached to the proton cf. Fig. 
7. The wavelengths after the dots are those of the series limits. 

TABLE II. Wavelengths in Angstroms of the Lyman, Balmer and Paschen of hydrogen with 
two Dirac monopoles attached to the proton cf. Fig. 8 and 6. 
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